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Abstract 

After considering the kinetical interaction principle (KIP) introduced in ref. 
Physica A 296, 405 (2001), we study in the Boltzmann picture, the evolution 
equation and the H-theorem for non extensive systems. The g-kinetics and 
the K-kinetics are studied in detail starting from the most general non linear 
Boltzmann equation compatible with the KIP. 
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I. INTRODUCTION 



It is common in equations of Boltzmann type to treat the collisions as local and instan- 
taneous. The spectral particle function in energy and momentum is given by the Dirac 
function d{e — / 2 m) . 

Assumptions to set Boltzmann equation are based on three main features: 

1. Triple coUisions are neglected, BBGKY hierarchy is truncated at the level of two- 
particle distribution function. One is limited to study diluite gases of neutral particles 
with short-range forces of interaction range much shorter than average particle sepa- 
ration and mean-free particle path between collisions. All the processes are governed 
by binary collisions; 

2. Two-particle distribution function is factorizable with one-particle distribution, i.e. 
two-body collisions are statistically independent. One does not retain information from 
previous encounters and memory about dynamical correlations; Boltzmann equation 
can be solved with a class of dynamical information. 

3. The particle kinetics underlying the Boltzmann equation is linear, i.e. the transition 
probability of a particle from one site to another is given by factors proportional to 
the distribution function / rather than to powers of / (as will be discussed in the case 
of non linear kinetics). 

A non local and non instantaneous treatment of binary collisions, as required in many dif- 
ferent physical problems related, for instance, to electron transport, Lorentz dense gas, rate 
properties in non ideal plasma and heavy-ion reactions, is strictly related to quasi-particle 
features and a Lorentz profile 5^(e, p) of the spectral particle distribution in energy and mo- 
mentum can be assumed. A complete and unique formulation and treatment of Boltzmann 
equation, in this case, is not yet available, particularly for systems far from equilibrium 
(which concept of quasi-particle works better in kinetic equations for non equilibrium sys- 
tems?) . 
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With the purpose of treating the Boltzmann equation in this context in the near future, in 
the present work we hmit ourselves to discuss a generahzation of the Boltzmann equation 
in order to treat non-extensive systems composed by particles whose kinetics is non linear 
and subject to binary, local and instantaneous collisions. 

This is achieved by employing the kinetic interaction principle (KIP) which is local and in- 
stantaneous but contains a generalized exclusion-inclusion principle which eventually takes 
into account of a non linear kinetics. The KIP has been introduced and described in refs. 
IHJ^. It permits to treat different statistical distributions already known in literature in a 
unified way and in time dependent conditions. 

The nonextensive statistics introduced by Tsallis has been considered in time dependent 
conditions in the framework of the Fokker-Planck pictures in refs. . Recently, the Tsallis 
kinetics has been studied also in the Boltzmann picture in refs. 

In the present effort, starting from the KIP and after introducing the nonlinear kinetics in 
the Boltzmann picture we consider in detail the Tsallis statistics and the ^-statistics p|,|2|JT0 



II. KINETICAL INTERACTION PRINCIPLE 

Let us consider a particle system in the Boltzmann picture where only point, binary 
collisions occur: (A + Ai ^ A' + A[). The most general non linear kinetics is described 
through the evolution equation: 



d d 1 dV(x) d 

' V- 



f{t,x,v) 



dt dx m dx dv 
= / ' 1 (Tv', n(t, X, v' V, v[ Vi) — 7r(t, x, v ^ v', Vi v[) , (1) 

The transition probabilities are defined according to the Kinetical Interaction Principle 
(KIP) [|,g by means of: 

n{t, x,v^ v', ^ v[) = Tit, X, V, v' , «i, v[) 7(/, /') 7(/i, f[) , (2) 

being / = f{t, x, v), f = f{t, x, v'), fi = f{t, x, Vi), f[ = f{t, x, v[) while the function 
7(/, /') assumes the form: 



7(/, /') = «(/)K/')c(/, /')• 



(3) 



The first factor a(/) is an arbitrary function of the particle population of the starting 
site and satisfies the condition a(0) = 0. The second factor b{f') is an arbitrary function 
of the arrival site particle population obeying the condition 6(0) = 1. The third factor 
/') = c(/', /) takes into account that the populations of the two sites, namely / and /', 
can eventually affect the transition, collectively and symmetrically. 
The KIP imposes the following form to the Boltzmann equation: 

f{t, X, v) 



d d 1 dV{x) d 
dt dx m dx dv 



= I d-v'rvrrv[T{t, X, V, v', V,, v[)c{f, f)c{h, f[] 
Jn 

X Hf) Kf) a{f[) h{h) - a(/) h{f) a{h) h{f[)] . 



After introducing the two functions: 



B{f, /i, /', f[) = c{f, n c(/i, f[) h{f) h{h) Kf) h{f[) 



and 



Kf) 



(4) 



(5) 



(6) 



one can write Eq. (HI) under the form: 



d d 1 dVix) d 

h V — — 

dt dx m dx dv 



X 



fit, X, v) 

I rv'rv,rv[Tit, X, V, v\ V,, v[)B{f, A, /', f[] 
{ exp [In Kif) + In /t(/0] - exp [In «:(/) + In «:(/i)] } . 



(7) 



We remark that this evolution equation describes a very large class of non linear general- 
ized kinetics which includes also the standard linear kinetics obtained after posing a(/) = /, 
bif) = 1 and c(/, /') = 1. 

We recall that in the frame of the present non linear kinetics the mean value of a given 
quantity A is defined through: 



<A>= f cTxcTvAf ; / (rx(rvf = l. 
Jn Jn 



(8) 
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III. GENERALIZED ENTROPY 



From Eq.(|^, in stationary conditions, we have: 



In K{fs) + In K{f,,) = In «:(//) + In 



(9) 



therefore the quantity In K,{fs) is the colhsional invariant of the system. Taking into account 
that during the colhsions the particle number, energy and momentum are conserved, we can 
express the colhsional invariant as [|I|: 

■1 



lnK{fs) = -f3 



mv + V{x) — n' 



(10) 



with f3 = 1/kgT. This last equation, which defines the stationary distribution of the system, 
can be written also under the form: 

■1 



fail's) = ^exp|-/3 



mv + V{x) — n 



(11) 



with /? /i = /? /i' + In Z and Z given trough / d"'v dP'x fs = 1. 
Let us consider the functional /C(t): 

ifs) 



IC{t) = -k^ I (Txrv I df\n- 
Jn J K 



(12) 



Of course, the stationary distribution can be obtained also by maximizing the functional 
Kit)-- 

51C{t) 







f = fs 



(13) 



In ref. it has been shown that, if the condition: 

dKif) 



df 



> 



(14) 



is satisfied and / obeys the Boltzmann equation (^, the quantity — /C(t) is a Lyapunov 
functional: 

dJC{t) 



it 

Kit) < /C(oo) 



(15) 
(16) 



In order to introduce the entropy S(t) of the system, we pose: 



/C(t) = S{t) -k^(5{E- /i') , 

and observe that the energy E of the system: 

E = I (Txd^v \-niv'^ + V{x] 
Jn 12 



(17) 



/ 



(18) 



is a conserved quantity: dE/dt = 0, as well as the particle number J-j^d'^x d'^v f = 1. By 
comparing the two expressions of the functional /C(t) given by Eq.s (p^ and we obtain 
that: 



S = -k^ j^ci"a;c/"i; J df In K{f) 



> 



The entropy S{t) obeys the H-theorem: 

dSjt) 
dt 

S{t) < S{oo) 

as can be verified immediately if we take into account Eq.s (|15D-(|1 



(19) 



(20) 
(21) 



IV. TSALLIS KINETICS 



Let us consider the kinetics defined by fixing «(/) through: 



\n[Z K{f)] = In^iZ f) 



(22) 



where lnq(a;) = [x^ — !)/(! — g) is the Tsallis logarithm and / is a normalized function: 
J^d^x d'^v f = 1. The evolution equation becomes: 



d d 1 dV(x) d ' 

h V — — 

dt dx m dx dv 

xT{t, X, V, v', v,,v[)B{f, f,, f, f[) 
x{{Z'Z[)-' exp [\n,{Z'f')+\n,{Z[f[) 
-{ZZ,)-' exp [ln,(Z/) + ln,(Zi/i)]} 
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f(t,x,v) = I d''v'd''vid''v[ 
Jn 



(23) 



where Z' ^ Z\ and Z'^ are the partition functions related to the distributions /, /', j\ and 
f[ respectively. 

Being dK{f)/df > for the system described by Eq. (|23|), the H-theorem is satisfied and its 
stationary distribution is given by: 



fs = ^ expg { -P 



- mv'^ + V{x) — fi 
2 



(24) 



where expg(x) = [1 + (1 — g) x]^/*-^ '^^ is the Tsallis exponential while the partition function 
isZ = St^<Fx<Fv exp^{-/? [mv"^ /2 + V{x) - ^]}. 

The above Eq. (^) is the Tsallis distribution and can be obtained also from the varia- 
tional principle, defined by means of Eqs. (|12D and (|1^), which now assumes the form: 



— / (TxcTv 
df Jn 



2 — q I — q 2 — g 2 



/ = o. 



From Eq. (pS]) it results clear that the entropy of the system is given by: 
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and can be written also as: 



<ln,(Z/)>4 



2-q 



2-q 



(25) 



(26) 



s = ^s',^_!\zf]+ 



2-q 2-q 
being the Tsallis entropy defined through: 



Jn ^ — Q 



(27) 



(28) 



V. THE K-KINETICS 

In ref. it has been proposed the following one parameter deformation of the logarithm 
function: 



2k 



(29) 
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which reduces to the standard logarithm as k — 0, obeys the scahng law In^^j (x" 
m In^^^j (x) and presents the following power law asymptotic behaviour: 

1 



x^o+ \2k.\ 
1 



X— > + oo 



\2k\ 



X 



k\ 



(30) 
(31) 



We consider the kinetics defined by fixing k(/) as follows: 



ln[Z«:(/)] = ln^,j(Z/) . 



The evolution equation (|^) becomes: 



(32) 



d d 1 dV(x) d 

h V — — 

dt dx m dx dv 



f{t,x,v) 



I rv'rv,rv[T{t, X, V, v', vr, v[) Bif, f,, f, f[) 
Jn 

X {{Z'Z[)-' exp[ln^^^(Z7') + ln,.>(^(/()] 
-iZZ^)-' exp[ln^,^(Z/) + ln^,j(Zi/i)]} , 



(33) 



and reduces to the standard Boltzmann equation if we pose k — and B{f, fi, f, f[) = 1. 
We observe that dn{f)/df > for V/t G -R so that the H-theorem still holds. 
The stationary distribution of Eq. (B3) is given by: 



-(3 [^mv'^ + V{x) - 



(34) 



with Z = j^dT'xdT'v exp^^^ {—[3 + V{x) — /i]}, being the ^-exponential defined as 

the inverse function of ^-logarithm: 



exp^^j(x) = (^Vl + x"^ + KX^ 



1/k 



(35) 



The ^-exponential reduces to the standard exponential as /t — > 0, obeys the scale law 
[exp^^^ (a;)]™' = exp^ ^ ^{mx) and shows a power law asymptotic behaviour 



exp^ (X) ~ 



(36) 
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It is easy to verify that the distribution 
the appropriate constraint of the entropy: 



can be obtained after maximization under 



2k Jn \1 + K 



1 - K 



f 



l-K 



(37) 



which reduces to the standard Shannon entropy Sq = —k^ J d^x d^v f \n{Zf ) as the deforma- 
tion parameter k — > 0. The entropy and the entropy S defined in Eq. (|l^) are connected 
through Sk = S — kgln Z. The variational equation ( [13|) reproducing the distribution (|34D 
can be written under the form: 







(38) 



VI. CONCLUSION 

In the present contribution, we have studied the kinetics of a non linear system in the 
frame of the Boltzmann picture. The main points are the following: 

The introduction of the KIP, given through a particular expression of the transition prob- 
ability (@), permits to describe the time evolution of a non linear system and imposes, by 
means of the functional K,{f) in Eq. (P), its steady state as stationary solution of the evolu- 
tion equation. 

The KIP imposes the entropy form of the non linear system. Its expression is given in 
Eq. ([T9|) , and, as shown in Eq. (|20|), obeys to the H-theorem whenever the condition 



dn{f)/df > is satisfied. 

Finally, we have discussed within the formalism here developed, as a working example, the 
well known Tsallis statistics and the K-statistics, already introduced and described by one 
of us in ref. [Q. 

The evolution equation ( pSf ) is very close to the equation recently proposed by Lima, Silva, 
and Plastino within a kinetic foundations of Tsallis thermostatistics. In some points we 
differ from them. These authors consider a hard sphere particle gas and a g-coUisional term 
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written in terms of a difference of two g-exponent correlation functions (before and after 
collision). Our coUisional terms is a function deduced from KIP and is a difference of two 
standard exponents. 

In ref. is reported that other possibilities also leading to Tsallis distribution can be 
obtained if exp^(x) is substituted by other positive increasing functions Fq{x) such that 
limq^i Fq{x) = exp(x). An attempt to unify the two derivations of Tsallis evolution equation 
is in progress. 
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